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ACCURATE  ANALYTIC  APPROXIMATIONS  AND  NUMERICAL 


SOLUTIONS  FOR  THE  STRUCTURE  OF  QUASI-STATIC 
LASER  DRIVEN  ABLATION 

I.  INTRODUCTION 

Previous  numerical  simulations  at  NRL  have  shown  that  density  and  pressure  profiles  for 
an  accelerating  D-T  shell  and  ablation  layer  can  be  temporally  superimposed  with  very  little 
scatter.  Because  of  the  presence  of  these  long  lived  "quasi-static"  profiles,  we  felt  that  it  was 
important  to  model  the  equilibrium  solutions  for  a  laser  ablation  layer  both  analytically  and 
numerically  and  to  develop  a  detailed  understanding  of  the  nature  of  this  flow.  Facility  with 
such  a  model  enables  us  to  calculate  the  quantitative  dependence  of  the  shell  thickness  and 
acceleration,  the  peak  density,  velocity,  and  temperature  of  the  shell,  the  width  of  the 
Rayleigh-Taylor  unstable  region,  and  the  distance  to  the  critical  surface  on  the  total  plasma 
mass,  the  critical  density,  and  the  absorbed  and  reflected  laser  flux.  This  "quasi-equilibrium",  as 
determined  by  the  input  parameters  of  the  configuration,  can  be  used  as  an  inexpensive  input 
or  driver  for  other  numerical  studies  and  as  an  equilibrium  solution  for  stability  analyses.  Fol¬ 
lowing  is  the  detailed  description  of  the  laser  ablation  layer  included  in  the  numerical  model. 

II.  STRUCTURE  OF  THE  ACCELERATING  SHELL 

AND  BLOWOFF  PLASMA 

For  our  analyses  we  use  the  basic  one-dimensional  steady  state  fluid  equations  for  a 
plasma  slab  transformed  into  an  accelerating  frame  of  reference 

continuity:  —  (pv)  -  0  (la) 

ax 
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The  isothermal  mach  number  is  defined  by 


(5) 


Only  two  of  the  three  equations  (2)  are  independent.  Depending  on  the  circumstances, 
we  found  that  various  pairs  of  the  three  equations  were  more  convenient  to  analyze  than  oth¬ 
ers.  Under  this  transformation  the  problem  was  reduced  to  solving  for  two  dependent,  dimen¬ 
sionless  variables  0  (the  temperature)  and  t)  (the  density)  based  on  one  independent  spatial 
variable  £  and  two  dimensionless  coefficients  N0  (primarily  the  thermal  conductivity  usually  a 
large  number),  and  M0  (the  isothermal  mach  number  at  the  density  peak,  usually  a  small 
number).  We  utilized  this  dimensionless  set  of  equations  (2)  to  generate  either  general  numer¬ 
ical  or  analytic  solutions  and  then  matched  these  solutions  to  the  specific  physical  constraints 
and  boundary  conditions  to  determine  the  desired  self-consistent  profiles.  The  numerical  solu¬ 
tions  were  obtained  by  a  Runge-Kutta  integration  of  the  non-linear  coupled  set  of  differential 
equations  (2)  while  the  analytic  solutions  were  obtained  by  making  appropriate  approximations 
to  the  equations  in  the  various  regions  and  then  solving  the  simplified  set. 

For  the  purpose  of  obtaining  general  analytic  solutions,  the  set  of  equations  (2)  is  con¬ 
veniently  divided  into  three  regions.  The  first  region  is  the  accelerating  shell  (£  <  0).  the 
second  Tegion  is  the  vicinity  of  the  density  peak  (£  =  0),  and  the  third  region  is  the  blowoff 
plasma  (£  >0).  Because  of  certain  basic  physical  trends  for  the  solutions  p  and  rj  in  each  of 
these  regions,  it  is  possible  to  make  simplifying  assumptions  that  allow  analytic  solutions  to  the 
equations  (2).  As  will  be  discussed  later,  there  are  other  admissible  solutions  to  these  equa¬ 
tions  with  counterintuitive  physical  trends  but  as  yet  these  have  not  been  sufficiently  investi¬ 
gated  to  determine  whether  they,  too  are  physically  real. 
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a.  The  Region  of  the  Accelerating  Shell  ({  <  0) 

Rewriting  equation  (2a)  into  a  more  convenient  form  we  find  that: 

5(0-1)  +  M0Hv~l)  -2{-  jj-  ^5/2-^-1|-  (6) 

The  traditional  physical  profiles  for  this  region  require  that  <  0,  0  <  1,  and  r)M„  ^  0. 

In  words,  the  density  and  temperature  should  decrease  away  from  the  density  maximum  at 
{  -  0  toward  the  inner  regions  of  the  accelerating  shell  (£  <  0),  and  for  this  region  the  flow  is 

isothermally  subsonic.  Therefore  from  (2b)  ^  1  and  correspondingly 

(7) 

di 

From  (6)  and  (7)  we  then  obtain  the  following  two  inequalities 

5(0-1)  +  M2(t,-1)  -  2£  <  0  (8a) 

and 

5(0  - 1)  +  M02(V  ~  1)  -  2{  +  ~  >  0.  (8b) 

Except  for  a  very  narrow  region  where  |( |  =  (N0  is  usually  large),  Eqs.  (8a)  and  (8b) 

imply  that 

5(0-1)  +A//(t,-1)-2£  =  0  (9) 

is  quite  a  good  approximation.  By  differentiating  (2a)  we  find  that 

5i£  +  M7  4jl-.2-  —  —  L5'2  -^-1.  no) 

di  0  di  N0  di  di  \ 
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Again,  except  where  |{|  is  on  the  order  of  ~  and  0S/2  is  varying  rapidly,  we  have  the 
approximation  that 


(11) 


Combining  Eqs.  (2b)  and  (11)  gives  us  the  further  approximation  that 


(tj/33)  *  0 


(12a) 


or 


t)/83  =  1.  (12b) 

This  latter  result  can  also  be  obtained  from  equation  (2c). 

Collecting  our  results,  we  have  as  a  general  approximate  solution  to  the  set  of  equations 
(2)  for  the  region  {  <  0, 


5(0  —  1)  +  M?(.y  —  l)  —  2{  ~  0  (13a) 

t)03  =  1.  (13b) 

It  is  interesting  to  note  that  this  approximation  is  independent  of  N0  and  therefore  the  flow  and 
profiles  in  the  accelerating  shell  are  not  basically  governed  by  thermal  conductivity.  Also  equa¬ 
tion  (13b)  states  that  the  entropy  P/py  is  approximately  constant  throughout  the  accelerating 
shell  which  behaves  like  an  adiabatic  gas. 


From  equation  (13b)  we  obtain 


dy  3_  dji 

dl  04  d{‘ 


(14) 
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Combining  this  result  with  equation  (11)  yields 


d£  * _ _ 

di  5/34  -  3  M} 


It  is  important  to  note  that  equation  (IS)  does  not  hold  when  {  -  0  and  -  1  |{  -  0  is  the 

di  1 

turning  point  for  the  density.  There  -  0  and  from  (2b)  -  ll.  Therefore  it  is 

dC  c-o  <*(  t-o  1 


evident  that  varies  from  1  to  approximately 


5  -  3A/„J 


very  quickly  over  a  narrow  region. 


It  is  possible  for  the  flow  to  approach  an  isothermal  sonic  point  (M„  r\  -  0)  in  this  region. 
From  (13a)  and  (13b)  this  occurs  when 

/3({,)  =  VK  (16a) 

and 

_ _  5  +  A/„2  ,  . 

{,  =  ly]W0 - Y^.  (16b) 

Substituting  (16a)  into  (IS),  we  have 


=  j 

di  l-c, 


where  (s  is  the  location  of  the  isothermal  sonic  point.  This  is  an  important  result.  It  implies 
that  equation  (2b)  may  not  be  singular  at  the  isothermal  sonic  point.  In  fact,  from  (14)  we 
expect  that 


dr)  __  _  3 
|c -t.  Mo 

Applying  L’Hospital’s  rule  to  the  full  set  of  equations  (2)  in  the  limit  where  /9  — *  M„i)  and 
— -•  1  yields: 
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lim 

H-Mh  di 
</{ 


Nn 


P 


3/2 


-  1 


2  Ml 


1  - 


4^-5 

q}/2 

1-- - IT 


P 


in 


-  \ 


1/2 


(18) 


N. 


Since  >>1  generally,  Equation  (18)  tells  us  that 


djL 

di 


3 

M} 


(19) 


as  expected.  The  fact  that  the  flow  approaches  a  well  behaved,  continuous  isothermal  sonic 
point  is  very  helpful  since  it  gives  us  a  convenient  position  to  define  as  the  inner  edge  of  the 
accelerating  shell.  This  definition  is  reasonable  since  no  physical  information  about  the  flow 
could  be  transmitted  across  this  point. 


b.  The  Vicinity  of  the  Density  Peak  (£  =  0) 


The  simplest  approach  to  determining  an  analytic  approximation  in  this  region  is  to 
expand  the  system  of  equations  (2)  in  a  Taylor  series  about  £  «*  0.  Doing  this  for  small  £  yields 
to  third  order 


0~  l+{  + 


-£(3tf0-5)  +  -£ 


NA3N0  —  5)  5—3  M02 


1  -  Mo 


-5(3A/0-5)  -  -y- 


(20a) 


and 


V 


1  + 


0/2 

\-M} 


(3 N0  -  5) 


£76  iV0(3iV0  -5) 
1-Atf  2 


5-3  Ml 

l— Mg 


ON0-5) 


7-5  M2 
l  -  M2 
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(20b) 
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Since  £  -  0  is  to  be  a  density  maximum 

>  0  (21) 

d?  Uo 

which  by  differentiating  equation  (20b)  twice  yields 

\0  >  j-  (22) 

This  is  not  a  very  stringent  criteria  for  N0.  Generally  for  physical  problems  of  interest 
N0  >>  5/3.  It  is  also  important  to  note  that  the  region  where  the  Taylor  expansions  (20)  are 
accurate  is  exceedingly  narrow.  Trying  to  utilize  these  results  for  even  a  distance  on  the  order 
of  the  thickness  of  the  ablation  layer  does  not  yield  valid  results.  In  fact,  equations  (2)  are 
included  more  for  completeness  than  for  their  practical  importance. 

c.  The  Region  of  the  Blowoff  Plasma  (£  >  0). 

Two  very  important  regions  lie  in  the  blowoff.  The  first  is  the  Rayleigh-Taylor  unstable 
region  where  the  gradients  of  density  and  pressure  have  opposite  sign.  The  second  is  the  vicin¬ 
ity  of  the  critical  surface  where  the  laser  deposition  occurs.  In  this  region  is  increasing 

di 

rapidly  from  its  value  of  1  at  £  -  0.  Corresponding,  /3  >  >  1  except  very  near  the  density 
peak.  It  is  expected  that  this  region  be  isothermally  subsonic  at  least  until  the  critical  surface  is 
reached.  Therefore, 

|M02(t,-1)  -2£|  <<  5(0-1)  (23) 

is  a  valid  approximation  except  in  a  very  narrow  region  about  £  -  0  where  it  is  in  error  by 
about  40%.  Here  retaining  £  is  important  and  the  approximation  is  more  properly  written  as 

A/02(t,-1)  <<  5(0-1)  -2£.  (24) 
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The  error  in  (23)  potentially  could  be  as  great  as  20%  in  the  vicinity  of  the  critical  surface  as 
well.  In  practice  (numerical  integrations)  it  never  becomes  that  large.  Under  these  approxima¬ 
tions  equation  (2a)  becomes  either 


di 


__  No  5 

2 

2  „ 

2  p^ 

0  - 

5  N0 

-  5  « 

(25a) 


or 


dp  __  N0  s 
di  2 


(25b) 


Equation  (25b)  can  be  solved  analytically  to  yield: 


£ 


\(PS/1-1)  + 


(pin  —  1)  + 


(pU2-l) 


Equation  (25a)  can  also  be  solved  analytically  to  yield  an  approximation  near  £  —  0  that  is  more 
exact  than  (26),  but  little  is  gained  by  including  it  here  due  to  its  complexity.  Though  it  is  not 
crucial  to  use  this  more  exact  approximation  to  define  the  Rayleigh-Taylor  unstable  region,  its 
use  does  improve  the  results  in  a  measurable  way.  For  £  >  0  it  is  advantageous  to  write  equa¬ 
tion  (2b)  in  terms  of  the  pressure  where 


-L.J- 

Po  Vt) 


(27) 


9 


J  II  GRI  NS 


10 


NRL  MEMORANDUM  REPORT  4167 


From  (28b)  we  obtain  two  interesting  results.  First 


and  second 


0  when 


dr\  2 

di  “  M}  ' 


Combining  (32)  and  (2b)  yields 


or  from  (5) 


0  when 


di 


M02v 


0  when  A/2-^  =>  1. 
dt 


(31) 


(32) 


(33) 


(34) 


Equations  (32)  and  (33)  or  (34)  give  us  respectively  the  rate  of  change  of  the  density  and  the 
temperature  at  the  turning  point  of  the  pressure.  It  is  just  this  offset  of  the  density  and  pres¬ 
sure  peaks  that  bounds  the  Rayleigh-Taylor  unstable  region,  so-called  because  density  and 
pressure  gradients  are  opposed. 


Combining  (33)  with  equation  (2a)  yields  an  equation  defining  the  relationship  of  the 
variables  at  the  pressure  peak 

_ N0  1 

Mlv  "  2  f}ifl 

The  relevant  question  is  whether  equation  (35)  has  a  solution  for  all  values  of  M0  and  N„. 
From  Eq.  (32),  in  the  Rayleigh-Taylor  unstable  region 

06, 


5(/3  —  1)  +  +M0H  TJ-D-2C 


(35) 
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The  only  other  important  modeling  necessary  for  the  region  {  >  0  is  near  the  critical  sur¬ 
face.  Here  we  integrate  equations  (1)  across  the  critical  surface  to  obtain  the  jump  conditions 
and  utilize  the  density  shelf  criteria  of  Felber'  and  Lee  at  a!.1  where 
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\AT  “  lyfih 


(i -My 

Ml  -  In  (Ml)  -  1 


(43) 


with  the  subscript  minus  referring  to  the  values  just  interior  to  the  jump  and 

i 


Vc 


—  (p,  is  the  critical  density). 

Pc 


(44) 


We  force  the  temperature  to  remain  constant  across  the  jump  and  the  heat  flux  to  be  zero  out¬ 
side  the  critical  surface.  We  did  not  believe  that  it  was  important  to  model  the  underdense  gas 
outside  the  critical  surface  since  the  density  is  so  low  there.  Felber1  has  chosen  to  model  this 
region  as  an  adiabatic  gas  and  that  would  probably  be  as  good  a  choice  as  any  other.  We  do 
require  that  the  flow  out  of  the  critical  surface  be  supersonic  and  that  gives  us  an  additional 
constraint  on  M„  that 

M,j  <  Ml  <  y  (45) 


This  constraint  is  also  included  in  Figure  1.  When  M„  and  N„  are  chosen  within  the  bounds 
depicted  in  Figure  1  the  solutions  to  equations  (2)  are  profiles  that  can  be  given  physical  explai- 
nations.  For  values  outside  that  regime  equations  (2)  are  still  solvable,  but  it  is  not  clear 
whether  the  profiles  generated  also  have  a  physical  significance  or  are  merely  mathematical 
curiosities.  This  question  is  worthy  of  further  investigation.  There  are  two  possible  approaches 
to  obtain  profile  solutions  to  equations  (2)  for  given  specific  boundary  conditions.  First  the 
coupled  set  of  non-linear  first  order  differential  equations  can  be  integrated  numerically  by  a 
Runge-Kutta  method  and  then  the  boundary  conditions  matched.  Second,  the  general  analytic 
approximations  (13),  (26),  and  (30)  can  be  utilized  instead  of  a  numerical  integration  and  then 
matched  to  the  specific  boundary  condition.  Conceptually  there  is  not  much  difference  in  either 
method.  Because  of  the  complexity  of  the  analytic  approximations  both  methods  require  the 
use  of  a  computer.  The  major  advantage  of  the  analytic  approximation  approach  is  that,  by 


13 


J  H  ORtNS 


allowing  an  error  in  the  results  of  a  few  percent,  significant  computational  time  can  be  saved. 
Also  the  analytic  approach  often  gives  a  better  insight  into  what  is  physically  occurring.  The 
three  available  boundary  conditions  not  previously  taken  into  account  are  the  total  plasma  mass, 
the  absorbed  laser  flux,  and  the  reflected  laser  flux.  Certainly  the  original  slab  mass  must  be 
encompassed  by  the  density  profile.  For  our  model  we  assume  that  all  the  mass  lies  between 
the  interior  isothermal  sonic  point  and  the  critical  surface.  Because  the  densities  outside  this 
region  are  exceedingly  small  compared  to  the  peak  density,  this  approximation  should  have  lit¬ 
tle  effect  on  the  results.  Therefore,  if  m  is  the  total  slab  mass  per  unit  area, 


or  from  (4) 


m 


(46) 


(47) 


Evaluation  of  the  integral  in  (46)  or  (47)  follows  from  an  integration  of  equation  (28b) 


or 


x: 


~^=-  +  M„  >/ij 

VTJ 


«c 

t. 


and  since  Mq  ys  -  /3S 


-7=-  +  Mo  yR  ~ 

y/Vv 


-7=-  +  Mq  yR 
y/Vs 


s;:  -  -fjLr  +  Mo1  yR  -  2 Mi  yR 

or 

dC  -  +  Mo1  yR  -  2M0  yR  . 


(48a) 


(48b) 


(49a) 


(49b) 


14 


NRL  MEMORANDUM  REPORT  4167 


Applying  the  approximations  (16) 

•C  *  7^  +  ^ -  1Mo'A  (50) 

Of  course,  in  (50)  tj,  and  fic  are,  as  yet,  an  undetermined  parameters. 


At  the  critical  surface  we  handle  the  jump  conditions  and  laser  deposition  much  as  Felber1 
did  where 


4>sl+2Azii_1+2.  '.-4« 


P- 


p0M0rlnficMl 


(51) 


and  we  have  used  the  constancy  of  pv,  and  equation  (5).  It  is  also  true  from  reference  2  that 

-  (1  -  Ml  0)  (1  -  -  j)  (52) 

CpoToPc  <t> 


By  use  of  equations  (43),  (47),  (50),  (51),  and  (52),  we  are  able  to  determine  the  position  of 
the  critical  surface  for  a  given  M0  and  N0.  This  is  accomplished  by  beginning  at  £  -  0  and 
moving  in  the  direction  of  increasing  {  using  a  Runge-Kutta  solution  of  equations  (2)  or  using 
the  analytic  approximations.  At  each  point  we  assume  that  we  have  reached  the  critical  surface. 
Then  from  (47)  and  (50)  we  can  obtain  the  value  r0  such  that  all  the  plasma  mass  would  be 
contained  within  the  density  profile  to  that  point.  With  the  r„  equations  (51)  and  (52)  will 
determine  the  corresponding  p„.  This  value  of  p0  may  then  be  used  to  see  if  equation  (43)  is 
satisfied.  When  (43)  is  satisfied,  we  have  reached  the  critical  surface.  Also  we  have  the  values 
of  p0  and  t0  needed  to  transform  the  general  profiles  into  the  specific  profiles  for  the  given 
boundary  conditions.  The  other  values  of  interest  v0  and  gcan  be  found  from  (4). 


The  preceeding  discussion  has  masked  a  rather  curious  fact.  Even  with  all  the  criteria  we 
have  specified,  we  cannot  determine  a  unique  solution  based  on  a  single  pair  of  M0  and  N0  but 
rather  have  obtained  a  whole  family  of  possible  "quasi-equilibrium”  solutions  based  on  values  of 
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M„  and  N0  only  constrained  as  in  Figure  1.  Maybe  we  have  not  specified  all  the  available 
boundary  conditions  but  it  is  difficult  to  think  of  another  two,  only  based  on  the  equilibrium, 
that  we  are  at  liberty  to  specify.  More  likely,  what  the  results  are  telling  us  is  that  there  really 
are  various  possible  "quasi-equilibrium"  states  and  just  which  one  would  occur  in  a  given  physi¬ 
cal  situation  might  be  determined  by  the  neccessarily  time-dependent  evolution  leading  up  to 
that  equilibrium.  Such  an  effect  might  be  an  initial  shock  which  impacts  an  entropy  to  the 
shell.  Others  may  exist. 

III.  NUMERICAL  RESULTS 

Figure  2  displays  the  profiles  of  density,  temperature,  pressure,  and  isothermal  mach 
number  for  the  "quasi-equilibrium"  where  Km, 7/2  -  10-33,  pc  -  4xl(T3  */cm3,  la  -  10.2  TW 
/cm2,  l,  -  10.2  7’W7cm 2,  m  —  .ITSmg/cm2  M0  =  .07,  and  Na  ~  700.  These  profiles  were 
obtained  by  a  Runge-Kutla  integration  of  the  equations  (2),  but  curves  generated  through  the 
analytic  approximations  discussed  earlier  are  virtually  indistinguishable.  These  are  rather  typical 
results  and  while  variation  of  M0  and  N0,  within  the  constraints  of  Figure  1,  does  change  the 
profiles  quantitatively  there  is  not  much  qualitative  variation.  The  density,  temperature,  and 
pressure  are  seen  to  decrease  smoothly  away  from  the  density  peak  toward  the  isothermal  sonic 
point.  For  this  case  the  thickness  of  the  shell  is  about  35  microns  and  the  peak  density  is  about 
.3g/cm3.  The  corresponding  temperatures  and  pressures  at  the  density  peak  are  respectively 
approximately  8  x  1012  cm2/sec2  and  2.5  x  1012  dyne/cm2.  The  acceleration  for  the  shell  is 
about  4  x  1015  cm/sec2.  It  is  impossible  to  see  in  the  figure,  but  the  offset  of  the  pressure  and 
density  peaks  is  about  .1  micron.  The  peak  pressure  is  only  slightly  greater  than  the  pressure  at 
the  density  peak.  In  the  blowoff  plasma  the  density  drops  sharply  while  due  to  the  thermal  con¬ 
ductivity  the  temperature  increases  rapidly.  It  is  interesting  to  note  that  the  decrease  in  density 
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is  nearly  balanced  by  the  increase  in  temperature  such  that  the  pressure  remains  relatively  con¬ 
stant.  The  critical  surface  is  approximately  100  microns  from  the  density  peak  and  the  critical 
temperature  is  about  4.5  x  1014  cm2/sec2.  The  density  and  isothermal  mach  number  just  inside 
the  critical  surface  are  respectively  approximately  4.5  x  103  g/cm3  and  .6.  These  profiles  have 
been  used  as  initial  conditions  for  time  varying  models  and  are  seen  to  remain  relatively 
unchanged  for  long  periods  of  time.  Therefore  they  are  truly  "quasi-equilibrium"  solutions. 
Further,  since  the  density  and  pressure  gradients  in  the  constant  adiabat  shell  are  aligned,  the 
shell  interior  is  Rayleigh-Taylor  stable. 

IV.  CONCLUSIONS 

Development  of  this  "steady-state"  model  is  still  continuing.  The  most  important  question 
yet  to  be  resolved  is  how  to  choose  M„  and  N„.  What  is  desired  is  to  relate  M„  and  N„  to  some 
physical  properties  that  would  remain  constant  for  the  solutions  while  external  parameters  like 
the  plasma  mass  or  the  laser  intensity  are  slowly  varied.  One  such  possibility  is  the  entropy  of 
the  accelerating  shell.  We  have  utilized  one  of  our  "quasi-equilibrium"  solutions  as  an  initial 
condition  to  a  time  dependent  simulation  model  and  followed  the  burning  of  the  shell.  During 
the  evolution,  the  shell  entropy  did  remain  relatively  constant.  We  have  also  examined  the 
constancy  of  the  shell  momentum,  but,  as  yet,  have  not  reached  any  definite  conclusions. 

Our  model  has  been  used  as  an  initializer  to  time  dependent  ablation  models  and  as  a 
pressure  driver  for  a  laser-target  compression  simulation.  For  this  latter  case  the  "steady-state" 
model  was  used  to  relate  the  variation  of  the  laser  intensity  to  the  pressure  at  the  surface  of  the 
shell.  Generally  we  used  the  analytic  approximations  for  these  studies,  since  it  entails  a  vast 
savings  in  computational  time.  We  also  plan  to  utilize  our  model  as  an  equilibrium  for  a 
piggy-back  analysis  of  Rayleigh-Taylor  stability  in  the  ablation  layer  plasma. 
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is  no  turning  point  in  the  pressure  profile.  For  the  region  to  the  right  of  the  curve  the  flow  out  of  the  critical  surface  is 


not  supersonic. 


•■■at  TV  TCNMaATVM 


Typical  solution  profiles.  Results  are  for  the  case  where  Km,n  -  10"3\  pr  -  4  x  1( 
TW/cm2,  /,  -  10.2  TW/cm2,  m  -  .773  mg/cm2.  Mn  ==  07,  and  =  700. 


